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Figure 1: Left: A/ ^ as a function of temperature for HISQ/tree, asqtad, and p4 at different Ni at light quark 
mass nil — O.OSmj. We compare with the stout action [|^]. Right: Chiral condensate M for the asqtad action 
parameterized by the scaling function fg. 



1. Introduction 

Understanding the properties of QCD at finite temperature is important for heavy ion collisions 
and the early universe. Since for the most part we need information about thermal QCD in the non- 
perturbative regime, lattice calculations are required. In the past couple of years we have developed 
a clearer understanding of the importance of cutoff effects in some thermodynamic quantities, and, 
as a consequence, the need for lattice calculations closer to the continuum limit and with improved 
actions with intrinsically smaller cutoff effects ^]. Here we report on the status of ongoing 
calculations of the HotQCD collaboration with two light and one strange quark using the rooted 
staggered fermion formalism. In particular we present new results for the asqtad action at smaller 
lattice spacing, i.e. with temporal lattice extent A^t = 12, and new results for HISQ ^ fermions 
with a tree-level improved gauge action (HISQ/tree) [Q] at A^^- = 6 and 8, which has intrinsically 
smaller cutoff effects. The strange quark mass m., is set to its physical value. For the light quark 
mass mi we keep the ratio ^ fixed, resulting, approximately, in lines of constant physics over the 
temperature range of this study [|l|]. The scale is set by using ri = 0.3106(8)(18)(4)fm |j5|], for the 
calculations with HISQ/tree action we use ro = 0.469fm to set the scale. The companion contribu- 
tion focuses on cutoff effects. In this contribution we investigate the chiral aspects of the finite 
temperature transition in 2 + 1 -flavor QCD. 



2. Chiral aspects of tlie quark gluon plasma (QGP) 

The light quark chiral condensate is defined as a derivative of the QCD partition function 
with respect to the light quark mass, = y^j^- It order parameter of the QCD phase 

transition in the massless limit but needs a multiplicative renormalization. For finite quark mass, 
it also requires an additive renormalization. The combination A/ v of light and strange quark chiral 
condensates at finite and zero temperature T cancels the multiplicative renormalization constant 
and the additive renormalization constant linear in the quark mass, 

Ais = — ■ — : . (2.1) 

( W)/,o - f ( W),,o 
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Figure 2: Disconnected chiral susceptibility for the asqtad action at different quark masses. On the right 
side Xi.disc has been rescaled by a factor yjini/nis. 

As can be seen on the left side of Fig. [T], the chiral condensate A/ ^ drops rapidly in the transition re- 
gion. As discussed in we observe that cutoff effects, particularly from taste-symmetry breaking, 
tend to shift the transition region to higher temperatures. Those effects are less pronounced for the 
HISQ/tree action. The transition region at mi = O-OSm^ for the new A^^ = 12 asqtad data is around 
T K, 170 MeV at this light quark mass, and for the new A^^- = 8 HISQ/tree data it is less. Let us take 
a closer look at the transition region. In this region physics is captured by the singular part of the 
partition function. For sufficiently small (light) quark mass mi the order parameter is described by 
a universal scaling function fc which depends on the critical exponents of the 0(2) (0(4)) univer- 
sahty class, the external fields mi and through H = and the reduced temperature AJ = 
see Ref. for more details. 

On the right side of Fig. |l| we plot the multiplicatively renormalized chiral condensate Mb = 
{YY)i and its parameterization given by the scaling function plus additional, scaling violating 
terms which stem from the regular part of the partition function, 

Mb{T,mi,m,)=h^/^fGit/h^^l^^)+atATH + biH, (2.2) 

where at and bi are constants, ^ = ^ and t = ^ with scales ho and to which need to be determined 
together with Tc in the chiral limit [^. The constants 5 and jS are critical exponents. Fits based 
on this parameterization yield good agreement with the data in the transition region as well as for 
larger temperatures. Less agreement is found for larger quark masses and at lower temperatures. 
The latter observation may lead to the conclusion that the scaling window in the hadronic phase is 
significantly smaller compared to the QGP phase. 

For a given lattice spacing the scaling analysis is especially useful for defining the crossover 
temperature Tp and for extrapolating that temperature to its value at the physical light quark mass 
or at the chiral limit of zero light quark mass (Tc). The crossover can be identified with the peak 
in the isosinglet chiral susceptibility, i.e., the peak in the derivative of the light quark chiral con- 
densate with respect to the light quark mass. From the scaling analysis we know explicitly how 
the crossover temperature depends on the light quark mass as it approaches the critical point. We 
will explore this approach in more detail in a future publication. For now we investigate the critical 
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Figure 3: Left: Xi. 

(jisc at fixed nij — O.OSnis for different A^t witli HISQ/tree and asqtad action. Right: The 
connected chiral susceptibility. 



temperature Tc by studying peaks in the chiral susceptibility directly. The "full" or isosinglet chiral 
susceptibility Xi"'^^ in terms of the light one-flavor chiral condensate given by 

xi"'^ =2^^{W)l = Xmsc + Xlcon with (2.3) 
;t:i,disc = ^{((TrD,-i)')-(TrD,-i)2} and Xlcon = ^i:,{Di\x,0)DY\0,x)) , (2.4) 

where D/ is the light quark Dirac operator and V is the spatial volume. We plot the disconnected 
chiral susceptibihty ;ti,disc for different A^-;- and m/, see left plot in Fig. For all parameters the 
disconnected chiral susceptibility develops a clear peak structure. The height as well as the position 
of the peak shows a sizable quark mass dependence. On the left plot in Fig. || we illustrate the cutoff 
dependence of ;ti,disc- We find a mild dependence on the lattice spacing in the height and position 
of the peak. For the HISQ/tree action these cutoff effects are reduced compared with the asqtad 
action. The peak of Xi.disc at finite A^^ and m/ defines pseudocritical temperatures Tp, which we use 
to determine the critical temperature Tc in the chiral and continuum limit. 

Before we present this analysis we consider, briefly, the quark mass dependence of ;t^i,disc- For 
sufficiently small quark mass the chiral condensate in the vicinity of the phase transition can be 
understood in terms of the three-dimensional 0{N) model. Below the crossover temperature and 
for sufficiently small light quark mass the quark mass dependence of the chiral condensate and its 
susceptibilities is controlled by contributions from Goldstone modes. These modes contribute a 
term proportional to ^/mj to the chiral condensate which, in turn, leads to a 1/ y/mj divergence in 
the chiral susceptibility. These features are evident in Fig. ^ 

In order to extract a continuum extrapolated critical temperature we first determine the pseud- 
ocritical temperature for several A^^ and m/ by fitting the peak position of Xi4isc with different fitting 
ansatze where we allow for an asymmetric shape. The different ansatze allow for an estimate of 
systematic error in Tp. The peak location in ;ti,disc is controlled by the singular part of the parti- 
tion function. In the chiral limit, the disconnected chiral susceptibility ;ti,disc at the pseudocritical 

temperature exhibits critical behavior, ;ti,disc(7p) ~ >nf , with the critical exponent 5. To extrap- 
olate the pseudocritical temperature we use an ansatz motivated by 0{N) models in leading order. 
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Tpimi,N^) = T, + b 



m 



+ c 



(2.5) 




where we use the critical exponent d = 
^ 0.54 from O(A^) model. The last term 
in Eq. accounts for the ^(c?^) cutoff 
dependence of staggered fermions. This 
ansatz is supposed to work well for suffi- 
ciently small quark mass and lattice spac- 
ings. For the fitting we take into account six 
values of Tp obtained with the asqtad action 
at A^T = 6,8, 12 at different m/. As can be 
seen in Fig. 0, the ansatz works reasonably 
well and we obtain a stable fit when omitting = 4 data. The black curve in this plot gives the 
continuum extrapolated (A^^ °°) pseudocritical temperatures obtained from the (combined) fit. 
At physical mass parameters, — 77> we then obtain a continuum extrapolated pseudocritical 
temperature for the chiral transition. Our preliminary estimate for the transition temperature gives, 



0.05 0.1 0.15 0.2 0.25 

Figure 4: The pseudocritical temperature for six en- 
sembles at different A^T and m/. The colored curves cor- 
respond to the combined fit of the ansatz in Eq. 2.5. The 
black curve is the continuum extrapolation, A^^ — > °o. 



Tp = (164 ±6) MeV. 



(2.6) 



The error is a combined statistical and systematic estimate where also the error of the scale setting 
coming from ri is included. 

The connected chiral susceptibihty Xi,con is an integral over the scalar, flavor nonsinglet me- 
son correlation function, and thus is probing the thermal properties of the medium. Note that in 
the continuum limit Xi.con can diverge in the thermodynamic limit only if the 11^(1) symmetry is 
restored, which is not expected to happen at the QCD transition temperature. In fact, lattice cal- 
culations ^ indicate that the scalar screening masses develop a minimum at temperatures slightly 
above the transition temperature. Therefore, the connected chiral susceptibility is expected to show 
a maximum above the chiral transition temperature, even in the chiral limit. This is supported 
by our calculations by comparing the peak position of the disconnected to the connected chiral 
susceptibility, as shown in Fig. |3[ Note that ;^i,con at finite quark mass suffers from additive and 
multiplicative renormalization. This issue will be addressed in more detail in a future publication. 

3. Deconfinement and chiral symmetry restoration 

When discussing deconfinement one often investigates quantities like the light (1) and strange (s) 
quark number susceptibility, Xi,s = VTT^fjTY' energy density £, or the Polyakov loop. Beside 
the Polyakov loop, where a relation (if any) to the singular part of the partition function is unknown, 
the other quantities are also sensitive to critical behavior. However, criticality is less pronounced 
in the sense that these quantities do not diverge at Tc in the chiral limit, in contrast, e.g., to the 
disconnected chiral susceptibility. For example, in the chiral limit the temperature derivative of the 
quark number susceptibilities and the trace anomaly in the vicinity of are respectively given by. 



dT 



T-Tr 



and 



■3p)/T' 



T-Tr 



T-Tc 



(3.1) 
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Figure 5: Light (left) and strange (right) quark number fluctuations for asqtad and HISQ/tree action com- 
pared to the HRG model (black curve). 



Note that at the critical temperature T^- the slope is controlled by a contribution that arises from 
the regular part of the free energy, while its variation with temperature is given by the singular 
part. Since the critical exponent a is negative it appears to be very difficult to reliably extract 
information about Tc from Xi,s because contributions from the regular part of the partition function 
may hide critical behavior High statistics lattice data would be needed to discriminate the singular 
from the regular contributions. Therefore, critical behavior in state-of-the-art lattice calculations is 
best seen in quantities like the disconnected chiral susceptibility. 

Let us take a closer look at Xi and Xs- For low temperatures the relevant degrees of freedom 
which contribute to Xi,s are hadrons and Xi,s is small. At high temperatures the degrees of freedom 
are quarks and gluons. The quark number fluctuations are then related to the number of quarks and 
increase with temperature towards the ideal gas limit. In Fig. |5| we compare Xi to Xs for different 
A^T and actions. We observe that in the vicinity of the pseudocritical temperature Xi as well as Xs 
appear to be rather smooth functions, the change in Xi, however, is more pronounced. In terms of 
the above discussion this is related to the smaller slope parameter compared to c[. The relative 
magnitude of these regular terms can quite naturally be understood in terms of a hadron resonance 
gas (HRG) model. A more detailed investigation of this issue is currently ongoing. Note that 
our calculations of deconfinement and chiral quantities indicate that deconfinement and the chiral 
transition appear at about the same temperature, see Fig. ^ and Fig. ||. 

In Fig. 1^ we show the trace anomaly for different actions and quark masses. As can be seen in 
the left plot, at large temperatures {e — 3p)/T'^ decreases towards zero. This is expected because 
at high temperatures the conformal symmetry is gradually restored. In the transition region and 
above {e — 3p)/T^ has its largest values. With respect to critical behavior the trace anomaly and the 



quark number susceptibilities share the same functional dependence, see Eq. ^ . Critical behavior 
thus is a subdominant feature in the trace anomaly, just like in the quark number susceptibilities. 
Therefore, it may not be surprising that {e — 3p)/T^ obtains its largest value at a temperature 
above the pseudocritical temperature. At temperatures below the pseudocritical temperature the 
trace anomaly decreases and, eventually, follows the prediction of the HRG model [^]. In the plot 
on the right in Fig. ^ we show the low temperature region of {e — 3p)/T^ in more detail. Compared 
with our earlier calculations we observe that our new = 12 asqtad and Nt = S HISQ/tree data 
appear to be closer to the HRG curve. However, it is also evident that at the current stage of our 
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Figure 6: Left: Trace anomaly for different actions at light quark mass m; ~ O.OSm, and the parameterization 
(black curve) from Ref. Right: Low temperature range of (e — 3p)/T^ at m/ = O.OSnis for HISQ/tree 
and asqtad action, and at m/ = O.lmj for asqtad action at A^^ = 8. We compare to the HRG model where 
resonances up to 2.5 GeV are included . 

calculations the data at low temperatures is not yet conclusive and statistics for the new data sets 
have to be improved. 

Outlook 

In this contribution we have given a status report of ongoing calculations within the HotQCD 
collaboration focusing on chiral properties at finite temperature for QCD with 2 + 1 flavors. We 
supplemented our earlier efforts with new calculations with asqtad fermions at A^.i- = 12 and m/ = 
O.OSnis, and with HISQ/tree fermions at A^^ = 6, 8 at nij = O.OS/n^. With this new asqtad calculation 
we extracted a preliminary continuum extrapolated transition temperature at the physical point, 
~ (164 lb 6) MeV. An improved determination of using universal scaling behavior is under 
investigation. The new HISQ/tree data in addition to upcoming HISQ/tree A^^ = 12 are much closer 
to the continuum limit and will enable us to further improve the continuum extrapolation. We 
remark that from our calculations we observe that the chiral transition and deconfinement appear 
at about the same temperature. Relating critical behavior to observables which are dominantly 
sensitive to deconfinement, i.e. the sudden change of degrees of freedom, is more subtle. 
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